Transport properties of clean and disordered Josephson junction arrays 
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We investigate the influence of quantum fluctuations and weak disorder on the vortex dynamics in 
a two-dimensional superconducting Berezinskii-Kosterlitz-Thouless system. The temperature below 
which quantum fluctuations dominate the vortex creep is determined, and the transport in this 
quantum regime is described. The crossover from quantum to classical regime is discussed and the 
quantum correction to the classical current-voltage relation is determined. It is found that weak 
disorder can effectively reduce the critical current as compared to that in the clean system. 
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The physics of the Berezinskii-Kostcrlitz-Thouless 
(BKT) transition in thin superconducting films is re- 
emerging as one of the mainstreams in the current con- 
densed matter physics. The interest is motivated by re- 
cent advances in studies of layered high-temperature su- 
perconductors the discovery of the superconductiv- 
ity at the interface between the insulating oxides [1, |3| , 
new studies in thin superconducting films uncovering 
the role of the two-dimensional (2D) superconducting 
fiuctuations 0, @| , and the intense developments in the 
physics of the superconductor-insulator transition where 
the BKT transition may play a major role 

The predicted benchmark of the transition that serves 
to detect it experimentally is the change of the shape 
of the I-V characteristics, V cx from a — above 

the transition, T > Tbkt, to a = 2{Tbkt/T) at T < 
Tbkt [1, @| ■ However, experimental data on supercon- 
ducting films show appreciable deviations from the theo- 
retical predictions and are still inconclusive ■ Among 
the possible sources of the deviation from the classic pre- 
dictions, one can consider the finite size effect, [HI, 
and effects of disorder 13|, [l5|. Another important 



issue is the role of quantum effects which become crucial 
when the BKT transition occurs at low enough temper- 
atures. In this Brief Report we will analyze the role of 
quantum effects in the BKT transition paying a special 
attention to the intermediate region of the interplay be- 
tween thermal and quantum contributions. We will dis- 
cuss the effect of disorder-generated vortices on the BKT 
transition, neglecting quantum fluctuations, namely the 
effective reduction of the critical current as compared to 
that in clean samples. 

Model. We choose a disordered Josephson junction ar- 
ray as a convenient discrete model for the 2D disordered 
superconducting film [l6| . The Hamiltonian describing 
the system is: 



2J 



J cos{(pi 



V'J-A^j) (1) 



where 



-2ei6j,k- We ignore single electron tun- 
neling and other sources of dissipation. The only non- 



vanishing elements of the capacitance matrix are its 
diagonal elements Cjj — AC (no summation over the re- 
peated index) and Cy — —C for the nearest neighbors 
i.e., the capacitance to the ground is assumed neg- 
ligible as compared to the mutual capacitances of the 
superconducting islands. The second sum in ([T]) is taken 
over all nearest-neighbor pairs on a square lattice. Ran- 
dom phase shifts Aij result from the deviations of the 
flux in a distorted plaquette from an integer multiple of 
the flux quantum $o = hc/2e [l7| . 

In the clean classical case, i.e. for Aij — and in the 
limit C oo, the physics of the system can be most 
adequately described in terms of vortices that experience 
the superconducting BKT transition at the temperature 
Tbkt — t^J 12, where J denotes the renormalized cou- 
pling constant. It is convenient to decompose the phase 
at the site i, as Lpi — (^["' -I- where Iv) and {sw) 

stand for the vortex and the spin wave part, respectively. 
Then, the vortex Hamiltonian can be written as 



Jtt rninij In 



(2) 



where Ef. denotes the core energy of a vortex. The sums 
are taken over the sites of a dual lattice; is the 
vorticity of the zth vortex, and we assumed that rui = 
0, where ^ denotes the superconductor coherence length. 

Next we want to include quantum fluctuations. After 
going over to the path integral representation of the par- 
tition function and integrating out the charge degrees of 
freedom, the action of the Josephson junction array in 
the hmit = /2C <^ J assumes the form [Ig, ^8\ 



S 



dr 



(3) 



The vectors ri(r) are the world lines of the vortices and 

Clean case. We begin with the discussion of a clean 
case. If we apply an external transport current, it will 
exert the force f j on the vortices, where j is the 



2 



current density [l9|. This generates an additional term 
— ^jTOif • in ([2]). In order to describe the effect of 
vortices on the current-voltage relation quantitatively, we 
consider the effect of vortices crossing the system trans- 
versely to the transport current. This motion dissipates 
energy. The Bardeen-Stephen flux flow resistance 
gives for the current-voltage {V — /) relation 



(4) 



where is the vortex density and p„ is the normal state 
resistivity. The rate equation for the vortex density is 



^2 

dtUy = r ^ — nl. 



(5) 



Here F denotes the rate of generation of free vortices, 
while the second term on the rhs of ([5]) describes their 
recombination, Tree denotes the recombination parame- 
ter. The steady state value n^, = {Trec^y^^ /S, of the 
vortex density determines the current-voltage relation. 

In order to determine F, we consider the appearance 
of a vortex-antivortex pair and its subsequent separation 
via tunneling or thermal activation under the influence 
of the external force f . In the clean case this process is 
symmetric, i.e., the coordinates of the vortex ri and the 
antivortex r2 satisfy ri — — r2 = r with f • r = fr. The 
action of the vortex pair can be rewritten as 



dr [MidrrY + U{r)] 



(6) 



where U(r) = 27rJln 



2/r + 2Ec. The problem 



effectively reduces to a single particle motion through 
one-dimensional potential barrier U (r) [2l| . 
The rate F is given by [2^ 



/•OO 

F~ / dEr{E)e-''/^, 
Jo 



(7) 



where T{E) denotes the zero temperature tunneling rate 
of a particle in the potential U{r) having an energy E. 
For low temperatures and hence E smaller than the bar- 



rier height Uq = 2ttJ 
WKB approximation is 



in(^: 



1 



+ 2Ec, T{E) in the 



TiE) 



(8) 



where ra/b{E) satisfy U{ra/b) — E (see Fig. [T]). In the 
following different regimes will be considered. 
(i)At zero temperature the only contribution in Eq. ^ 
comes from E = Q. The generated voltage for small cur- 
rents (/C/ J < 1) is 



V - Fi/2 
5(0,0) 
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FIG. 1; Potential barrier for the separation of the vortex- 
antivortex pair. 

ci is a positive constant of the order of unity and 
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(10) 



is the action of the classical path of the particle in the 
potential —U{r) with the energy E and mass 2M. The 
result ([9]) is in an agreement with that of Ref. 23| where 
it is obtained using the different technique [2^. We find 
that the result ([9]) holds also at finite temperatures as 
long as 



1 



1 



C2 Vn2JM 



(11) 



where C2 is positive constant of the order of unity, 
(ii) At intermediate temperatures To < T < T* , where 



T* = 



2-K 



-U"{rc) 
2M 



M / 1 
27r V MJt: ' 



(12) 



the main contribution in Eq. ([7]) comes from the station- 
ary point Et- Therefore, V ~ el^p[-S{ET,T)/2h]. Et 
depends on the temperature and is implicitly given by 
the equation 



- = 2VM 
T 



rb(-Br) 



dr- 



1 



Ta{ET) \JU{r) - Et 



t{Ej 



(13) 



where t{E) can be interpreted as the period of the clas- 
sical motion of a particle with the mass 2M and energy 
i?, in the potential —U{r). Since t{E) is the mono- 
tonically decreasing function of E for small currents, 
Eq. has the unique solution Et for every T in a 
range Tq < T < T* . We come back to the discussion of 
the voltage characteristic in this regime below. 
(iii)At even higher temperatures T* < T < Tbkt, the 
decay rate is dominated hy E > Uq 13, 25| and the ther- 
mally activated breaking of vortex pairs dominates the 
dynamics. Then, the decay rate is given by the Arrhe- 
nius law F ~ exp[-Sciass/h] where Sdass = hUo/T. The 
voltage-current relation reads [1, [2^ 



5{T) ^I + ttJ/T. (14) 



Taking into account the presence of other vortices by 
replacing J ^ J, the coefficient assumes a universal value 
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FIG. 2: Dynamic phase diagram in current-temperature co- 
ordinates showing different types V{j,T) dependencies for 
T < Tbkt- The dashed and the sohd hues sketch To(j) and 
T*{j), respectively. In the domain T <To quantum tunnehng 
of vortices dominates the vortex dynamics, while at T > T* 
the voltage-current characteristics is determined by the ther- 
mally activated motion. In the shaded region the quantum 
correction to the classical result, given by Eq. (|17|l . applies. 



S{Tbkt) = 3. 

(iv) At T > Tbkt a finite density of free vortices appears 
in an equilibrium, and the system is characterized by a 
linear current- voltage relation for small enough currents. 

Next, we consider crossover from the quantum- (T < 
Tq) to the classical regime (T > T*) in more de- 
tail. Within the semiclassical approximation the de- 
cay rate is given, with the exponential accuracy, by 
r ~ exp[— 5'inin/S], where S'min is the action of the trajec- 
tory minimizing the Euclidean action of Eq. ([6|) . For tem- 
peratures below To the extremal action is S'min = '5'(0, 0), 
in the intermediate region (Tq < T < T*) the minimal 
action is S'min = S{Et,T), and in the high temperature 
regime the trajectory extremizing the action is time in- 
dependent, and therefore S'min = KUq/T. We find that 
Smin at T* has a continuous first derivative with respect 
to temperature, while the second derivative has a jump: 



dS(i;T,T) 



dSc 



dT 

A^S{Et.T) 



dT 



dT2 



< 



dT2 



(15) 



Following Ref. [27| we call this situation a second-order 
transition at the crossover point [2^. The result of 
Eqs. is a general property of a massive particle 

trapped in a metastable state formed by a potential U (r) , 
provided t{E) is a monotonously decreasing function of 
energy [29| . 

Generally, in the case of a second-order transition the 
trajectory extremizing the action can be written as 27 1 



r(T) 



/27rT 
cos I ^ nr 



(16) 



where the coefficients |a„| <C |ai| {n > 1) are small 
near the transition temperature T*. Substituting r(r) 



in Eq. ([6]), the action can be expanded in powers of a„, 
yielding an effective action S « UoH/T+aal+f3af, where 
the coefRcient a is negative below T* and vanishes at the 
transition temperature T* 27|. Then the coefficient Oi 



can be found from the minimization of the action S and 
the minimal action is S'min = Uoh/T—a^/ (4/3). Following 
Refs. 27[ , we determine the coefficients a and /3 and find 
a quantum correction to the classical result of Eq. 



A = 



(T2 - T*2)2 ^MJ' 



hf l + 2(l-4(T/T*)2)-i' 



(17) 



This result is valid near the transition, for temperatures 
approaching T* from below, see Fig. (5] We conclude that 
quantum effects significantly enhance the decay rate in 
comparison to the classical rate for the asymptotically 
small currents. It would be interesting to probe the re- 
sult of Eq. (jl7p in experiments. 

Disordered case. Next we include disorder into the con- 
sideration, in the limit C —^ oo. The phase shifts are as- 
sumed to be uncorrelated from bond to bond, and each is 
Gaussian distributed with the mean value and the vari- 
ance 



A.,, = 0, 



42 



(18) 



respectively. Then, an additional term '^imiV{ri) is 
generated in where V{ri) = 27rJ^^ Qj ln(|ri — 

Tjl/O- Qi = (l/27r) I]<p/ag> ^rc the frozen charges 
sitting on the dual lattice in the center of a plaquette 
whereas the sum is over a plaquette formed by 4 bonds. 
From (HH]) follows {V{r^) - V{rj))^ = AnaJ^ \n{\r^ - 

It was shown in Ref. [l3| , that the system in the classi- 
cal case, at T = undergoes a disorder driven transition 
from the 'ordered' BKT state to a disordered phase at 
the critical disorder strength CTc = tt/S. In the ordered 
BKT phase vortices appear, on average, only in a form 
of the bound pairs. Indeed, the energy of a vortex pair 
with the separation R and mi — —7x12 = 1 in a clean 
sample is given by 27rJln(i?/^). Since V{Yi) is Gaussian 
distributed, the typical energy gain is —IJ^J-na ln(i?/^) 
which is smaller by a factor ~ (ln(i?/^))~^/2 than the en- 
ergy cost of a pair. However, the maximum energy gain 
of a vortex dipole in a region of linear size L > Ris larger 
by a factor V2\nN than the typical energy gain, which 
arises from the N independent realizations of the vortex 
positions [lj|. The disorder potential, that one vortex- 
antivortex pair of size R + dR feels, is uncorrelated when 
the pair is translated over a distance larger than R [sot . 
Therefore, we introduce a lattice with a lattice constant 
R. Since also the correlations of disorder potential inside 
the cell give only subleading-order corrections [soj, we 
estimate N w {L/ Rf{R/if{2TTR/C)dR/t For dR w i?. 
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we get the free energy of the pair at T = 



E Rr! 27rJln 



1 



Ua In {LR/e ) 



(19) 



Thus, if i? « L, the total energy of the coresponding vor- 
tex pair becomes negative and free vortices are favored 
by disorder provided cr > CTc = tt/ 8, in an agreement with 



the renormahzation group result in Ref. [13[. Note that 
strictly speaking these vortices are "pseudo-free" since 
despite the fact that their attraction is overruled by dis- 
order, they remain pinned by the same disorder-induced 
forces. It follows from the above reasoning that even for 
CT < CTc some rare vortex pairs of the negative energy 
can appear. From (|19p we get that their maximal size 
is i?c ~ £,{L/C) ^"a/'-i J which reaches the size of the sys- 
tem for cr — > (Tc — 0, as expected. Typically there is a 
single dipole of the size Rc in the system. If we divide 

the system into subsystems, each part will contain a 

1 

dipole of the maximum size Rm ~ RcM 2a^/„-i ^ xhe 
density of dipoles of the size Rm is C^'^{Rm /C)'^^'^^'^'"'^''^ 
at T = 0, in agreement with Ref. 30]. 

We further determine the critical current. If the trans- 
port current is strong enough, it will depin vortices such 
that the dissipation sets in. A crude estimate for the 
critical depinning force at T = and a < ac is given by 



fc 



J 



(20) 



since smaller dipoles are depinned at larger forces. The 
influence of disorder on the voltage-current relation is left 
for further studies. 

Conclusion. We have investigated transport properties 
of Josephson junction arrays taking into account the in- 
fluence of quantum fluctuations on the unbinding of vor- 
tex pairs for J. At sufficiently low temperatures 
the quantum tunneling of vortices turns out to be more 
probable than their thermal activation. We have derived 
the V-I relation corresponding to the quantum creep of 
the BKT- vortices and found the range of temperatures, 
< T < To, where this law is applicable. We have 
determined the temperature T* above which the ther- 
mally activated breaking of vortex pairs dominates the 
vortex nucleation. We have discussed the region of in- 
termediate temperatures Tq < T < T* where a crossover 
from classical to quantum behavior occurs, and found the 
quantum correction to the classical result, see Eq. (fTT]) . 
The results are schematically summarized in Fig. [2] and 
can be straightforwardly extended to the quantum limit 
Ec ^ J, where the transport is mediated by the motion 
of charges dual to the superconducting vortices, via the 
standard dual transformation. Moreover, in the presence 
of positional disorder and for C — > oo, we have shown 
that additional vortices generated by the disorder con- 
tribute to transport, effectively reducing the critical cur- 
rent as compared to that in a clean system. 
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